Roughness distributions for l/f a signals 
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The probability density function (PDF) of the roughness, i.e., of the temporal variance, of l/f a 
noise signals is studied. Our starting point is the generalization of the model of Gaussian, time- 
periodic, 1// noise, discussed in our recent Letter M, to arbitrary power law. We investigate three 
main scaling regions (a < 1/2, 1/2 < a < 1, and 1 < a), distinguished by the scaling of the 
cumulants in terms of the microscopic scale and the total length of the period. Various analytical 
representations of the PDF allow for a precise numerical evaluation of the scaling function of the 
PDF for any a. A simulation of the periodic process makes it possible to study also non-periodic, 
thus experimentally more relevant, signals on relatively short intervals embedded in the full period. 
We find that for a < 1/2 the scaled PDF-s in both the periodic and the non-periodic cases are 
Gaussian, but for a > 1/2 they differ from the Gaussian and from each other. Both deviations 
increase with growing a. That conclusion, based on numerics, is reinforced by analytic results for 
a = 2 and a — > oo, in the latter limit the scaling function of the PDF being finite for periodic 
signals, but developing a singularity for the aperiodic ones. Finally, an overview is given for the 
scaling of cumulants of the roughness and the various scaling regions in arbitrary dimensions. We 
suggest that our theoretical and numerical results open a new perspective on the data analysis of 
1/ ' f a processes. 

PACS numbers: PACS numbers: 05.70.Ln, 64.60.Cn, 82.20.-w 



I. INTRODUCTION 

The power spectra of fluctuations scales with frequency 
as S(f) ~ 1//" in a large variety of physical, chemical 
and biological systems Q . This power law behavior 1/ f a 
often persists over several orders of magnitude with cut- 
offs present at both high and low frequencies, and with 
typical values of a in the range 0.8 < a < 4 |. In a 
somewhat loose terminology, all these systems are said 
to display 1/f noise although good quality data with 
a very close to 1 exist only for the voltage fluctuations 
when a current is flowing through a resistor || |J. Phe- 
nomena with a =/= 1, however, are abundant, examples 
being the white-dwarf light emission || , the flow of sand 
through hourglass Q , ionic current fluctuations in mem- 
brane channels R, number of daily trades in the stock 
market ||, water flows of rivers 0, the spike trains of 
nerve cells [^Oj, the occurrence of earthquakes 11 , the 
traffic flow on a highway |l^, |l3), the electric noise in 
carbon nanotubes |Q and in nanoparticle films |L5| , the 
interface fluctuations [ fl6| |, dissipation in turbulent sys- 
tems |L7j, and the list could be continued. 

A well understood example of l/f a type behavior 
is the dynamic scaling observed at equilibrium critical 
points where the power-law correlations in time are gen- 
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erated by the infinite-range correlations in space. Most 
of the examples listed above, however, are related to non- 
equilibrium phenomena and a similar level of understand- 
ing does not exist. There have been many attempts at 
identifying possible generic mechanisms leading to scale 
invariant fluctuations, a notable example being the con- 
cept of self organized criticality jl8], |l9). It is clear, how- 
ever, that not all the systems showing a 1 / f a noise fit into 
a single scheme and, perhaps, at this stage one should 
pursue a less ambitious aim of developing more detailed 
characterizations of non-equilibrium universality classes. 

In equilibrium systems, the static universality classes 
are determined by the (i) dimensionality, (ii) symmetry 
of the order parameter, and (iii) the range of the inter- 
actions. Further specifying the conservation laws and 
the coupling of the order parameter to conserved quanti- 
ties defines then the dynamical universality classes. The 
statics and dynamics are, however, intertwined in non- 
equilibrium systems and an exponent in the 1 //" behav- 
ior carries information about both. Thus measuring a 
single exponent a (a characteristic situation when mea- 
suring time-series of noise) does not determine the uni- 
versality class of the system even within the framework 
of an equilibrium- type theory. Having a single exponent, 
however, one can still go on and try to ascertain whether 
two systems belong or not to the same universality class. 
This can be done by using the data to measure and com- 
pare the scaling functions associated with the finite-size 
scaling of some global physical quantity such as, e.g., the 
order parameter in a critical system or the roughness of 
an interface. 

The remarkable features of scaling functions are that 
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they are obtained without any fitting procedure and, fur- 
thermore, they usually converge fast as the system size 
is increased. Thus one can build a picture gallery that 
can be effectively used to identify systems belonging to 
a given universality class. Indeed, such an approach has 
been useful in establishing connections among rather di- 
verse processes such as massively parallel algorithms [^0| , 
interface dynamics in the d = 2 Fisher-Kolmogorov equa- 
tion pHl, dissipation fluctuations in a turbulence experi- 
ment [2^] and the interface fluctuations in the d = 2 
Edwards- Wilkinson model (XY model) pi p). This 



scaling function approach has also helped to clear up 
some questions about the upper critical dimension of the 
Kardar-Parisi-Zhang equation |26) . 

In our attempts to expand the picture gallery, we have 
recently derived |lj the following result: The scaling func- 
tion for the roughness-distribution of a Gaussian periodic 
1// noise signal is one of the extreme value distributions, 
the Fisher-Tippet-Gumbel distribution j2?|, . This is a 
rather unexpected and interesting result and we feel that 
it is important to explore its generality and limitations. 
First, because ideas about extreme statistics playing a 
role in strongly correlated, scale invariant systems have 
been much discussed pi], |2j| |9[ ||(| |3l| , and this result 
may provide a foundation to those speculations. Sec- 
ond, because this is the first instance when one of the 
distributions associated with extreme statistics emerges 
naturally, and in a mathematically precise manner, for a 
quantity which is not of extremal character a priori. 

There are three elements underlying the above extreme 
statistics results: Gaussianity, 1// power spectrum, and 
periodicity. Thus the problem can be generalized by con- 
sidering effects of (I) non-gaussianity, (2) a generalized 
power spectra of the form 1/ f a , and (3) non-periodicity 
(the experimentally realistic situation). In this paper, we 
shall concentrate on points (2) and (3) leaving the more 
difficult problem of non-gaussianity to a later study. 

One expects that changing the exponent a will change 
the roughness distribution of a signal and, indeed, there 
are analytical results jl], |3^, |3^| for a = 1,2,4 which 
demonstrate this explicitly for periodic signals (d=l di- 
mensional systems). Similar results exist also for higher 
dimensions where the roughness of a d=l signal is re- 
placed by the roughness of a higher dimensional inter- 
face [^4[ [34], ^5). Thus the questions we address in 
connection with the a-dependence is not about its ex- 
istence but about its magnitude (observability). More 
specifically, we ask if there is an interval of a where the 
a-dependence is absent or negligible from experimental 
point of view? 

The other problem of our concern is the question of 
the effect of the boundary conditions (BC-s). When an- 
alyzing a signal and building a distribution function, one 
usually divides the signal into equal "time" intervals and 
measures the quantity of interest in every interval. Thus 
the signal within a given interval is not periodic (usually 
the total signal is not periodic either) and the question 
arises whether the experimentally used BC, i.e., when 



one studies a "window" within the signal, would affect 
the scaling functions or not. This shall be a central issue 
in this paper, so we introduce the abbreviations PBC and 
WBC, for periodic and "window" BC, respectively. 

The dependence of scaling functions on the bound- 
ary conditions has been discussed in equilibrium critical 
phenomena and it is known that such dependence exist 
p6| , p7[ . Thus our aim here again is not to show its pres- 
ence but to gauge its magnitude for various a values and 
investigate for which values of a this effect is absent or 
negligible. 

As we show below, the scaling functions can be derived 
analytically in case of PBC. Similar analytical treatment 
for WBC were achieved only for a — 2 and in the a — > 00 
limit, otherwise we had to resort to simulations in order 
to obtain the finite alpha results. Observing the scal- 
ing functions (displayed on the Figures below) one can 
deduce the following qualitative trends as a is decreased. 

(1) The PBC scaling functions change smoothly apart 
from a w 1 where the function is sensitive to small 
changes in a since a singularity develops at a = 1. Even 
this singular behavior can be smoothed out, however, by 
an appropriate change of the scaling variable. The scal- 
ing function approaches a Gaussian as a — > 1/2 and it 
remains a Gaussian for a < 1/2. 

(2) The WBC scaling functions display strong a de- 
pendence at large a. The dependence becomes weaker for 
a < 2 and, similarly to the PBC case, the a-dependence 
disappears entirely for a < 1/2 where the function be- 
comes a Gaussian. 

(3) Comparing the PBC and WBC scaling functions, 
one can observe that their difference is large for a > 4, it 
is easily noticeable in the range between 1 and 4, while 
it becomes harder to distinguish the functions for a < 1, 
and the functions become identical Gaussians for a < 
1/2. 

In order to demonstrate the above observations, we 
shall introduce a model of 1//" noise in Sec. II where the 
roughness distribution, the appropriate choice of scaling 
variables, and some general scaling properties will also be 
discussed. Sec. [II contains the analytical calculations for 
the case of PBC while Sec. IV is devoted to presenting 



both the simulation and analytical results for the WBC 
case. We close with discussing generalizations to higher 
dimensions in Sec. [v| 



II. PROBABILITY DISTRIBUTION OF 
PERIODIC 1// Q SIGNALS AND THEIR 
ROUGHNESS 

A. Model for periodic, perfectly l/f a , Gaussian 
noise 

Noisy signals are in principle fully characterized by 
their path probability density functional. However, the 
most often used characterization is by the power spec- 
trum, whose experimental recording is straightforward. 
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This leaves several properties unspecified, such as the dis- 
tribution of the phases, and whether the noise is Gaus- 
sian. Furthermore, the usual Fourier expansion of the 
signal in an interval implies periodicity outside the in- 
terval. While that expansion can be used well for data 
analysis, the physical implication for the outer part of the 
signal, which influences by long-time correlations the in- 
ner part, is obviously non-realistic. Accordingly, the role 
of initial and final conditions in time, a property that we 
shall call boundary condition (BC), are generically ne- 
glected. Whereas it is arguable that in a stationary pro- 
cess, long time after the turning on of the experimental 
device, the boundary conditions should not have much 
influence, when a model is constructed for 1// Q noise, 
the boundary conditions should be defined with care. 

We now describe the main properties of the model 
treated in this paper. It should be emphasized that we 
do not study the underlying microscopic mechanism that 
may lead to l/f a noise, this has been done on many 
different levels of simplification for a large number of 
systems, for a review see S. Rather, we construct a 
stochastic model, where the noisy signal is given by its 
path probability density functional, having the generic 
properties of observed 1// Q noises. In a sense it will be 
a minimal model tailored to exhibit a few properties we 
prescribe. 

Firstly, the noise we consider is periodic. While this 
is obviously not valid in a real experimental situation 
over the time interval of the measurement, as we have 
seen recently [Q, it yields theoretical predictions not 
very far from what is measured in a real signal. Fur- 
thermore, the periodic boundary condition (PBC) allows 
for the straightforward numerical simulation of the noise, 
thus making it possible to study numerically the statis- 
tics in cut-out time intervals that are much smaller than 
the entire period. Such a "window" boundary condition 
(WBC) is, we surmise, a more faithful representation of 
experimentally realistic BC-s. Indeed, while the effect 
of the outer signal onto the inner part is important due 
to long-time correlations, probably it depends little on 
whether the signal was simply turned on long time ago, 
or is periodic on a time scale much larger than the "win- 
dow" . Looking at it from another angle we can say that 
the PBC is equivalent to saying that the signal is best ex- 
panded in the usual Fourier basis, thus what remains is 
to specify the probability distribution of the Fourier coef- 
ficients. Note that the noisy signal in a "window" should 
not be expanded in such a basis in theory, because the 
noise is manifestly aperiodic. 

Secondly, we restrict ourself to Gaussian noise and as- 
sume that the Fourier modes are independent random 
variables. The reasons for this working hypothesis are 
that many experimental data concern power spectrum 
measurements of Gaussian noises and provide no infor- 
mation about the coupling between the modes. 

Thirdly, we assume that the phase of each complex 
Fourier component is random and uniformly distributed 
in [0, 27!"]. Hence the probability distribution will depend 



only on the modulus. 

Last but not least, we consider perfectly 1// Q noise. 
Here we understand that the variance of the independent 
Fourier amplitudes of a path is a pure power with expo- 
nent a. Furthermore, there are no upper and lower cutoff 
frequencies other than the natural ones determined by 
the observation time and the microscopic characteristic 
time unit, respectively. 

The above ingredients specify the path probability of 
the stochastic time signal, this will be our model, as dis- 
cussed in detail below. (Note that it is a generalization 
of the case a = 1 to arbitrary a.) 

The stochastic trajectories h(t) are periodic, h(t) = 
h(t + T), thus the trajectory can be expanded as 

N 

h(t)= ^ c n e^ mt ' T , c* n = c_ n . (1) 

n=-N 

Here we have used a cutoff N meaning that the time 
scale is not resolved below r = T/N. This is needed to 
deal with some singularities and may have microscopic 
physical origin. Whenever possible we shall, however, 
take iV — ► oo, and even if this were a singular limit, 
we shall wind up with scaling functions not containing 
the microscopic time unit, so the value of r may be left 
undetermined. 

The stochastic properties of a signal h(t) are fully char- 
acterized by our specifying the probability density func- 
tional 

P[h(t)] = e - s W*)l, (2) 

or, equivalently, the probability density function (PDF) 
for the Fourier coefficients c„ for n — ... N as 

P({c n }) = Ae- S « c "». (3) 

Our model for periodic, perfectly l/f a signals is defined 
by the action 

JV 

S({c n }) = 2aT 1 - a Y, na M 2 > ( 4 ) 

n=l 

JV 

A = Yl 2(Tn a T 1 - a n- 1 , (5) 

n=l 

where the probabilistic variables are the real and imag- 
inary parts of Cn-s. 1 We assume translation invariance 



1 In time-representation the action is of the form 

/-i\ Q r T ( d a /' 2 h\ 2 
*[*(*)]= a, (^) I , 

where the fractional derivative is understood on the Fourier rep- 
resentation so that it acts upon the phasefactor as d^e lat /dt@ = 
(ia)Pe' lat . Dimensional analysis on Eq. (hi) yields the scaling by 
T used in (GJ). 
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in h space, therefore the action must not depend on the 
constant part Cq and we set that to zero hereafter. Note 
furthermore that the above action means that different 
Fourier modes are uncorrelated. The coefficient a makes 
S dimensionless and can be understood as the recipro- 
cal noise intensity parameter. Since we are after scaling 
functions, the value of a will not matter. The above ac- 
tion was set up so that the power spectrum is 



'•n\ 2 ) OC — 

n a 



(6) 



thus the process is indeed of the l// Q -type. 



B. Roughness of a signal 

Our aim is to characterize the signals by some global 
properties. As the time average cq in Eq. ([j]) was set to 
zero, the simplest quantity worth considering is the mean 
square width, i.e., roughness of the signal. This has been 
studied for the 1/f noise (a = 1) in Ref. 0, the Wiener 
process (a — 2) in Ref. [p2[ , and for curvature-driven 
interfaces (a — 4) in Ref. (33| . Its general definition is 
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h{t)-h{t) =h 2 {t)~h{t) 



(7) 



where overbar me ans t ime average over the entire period 
T. Note that now h(t) = 0. In Fourier representation we 
get 



w 2 



(8) 



hence w 2 is in fact the integrated power spectrum. 

The roughness w 2 is associated with a given h(t) tra- 
jectory, it varies when different instances of the paths 
are taken from their ensemble. Thus w 2 is a probabilistic 
variable, whose PDF can be expressed as 



P (w 2 ) = J S[w 2 - h 2 {t) + h(t) J P[h{t)}Dh{t). (9) 

In the more practical Fourier representation P (w 2 ) as- 
sumes the form 

P{W 2 ) = J S(w 2 - 2J2 P({Cn}) 
N 

x J^J dRec n dlmc„ (10) 



It is useful to introduce the generating function of the 
moments of P (w 2 ) as 



G(s) = / dw 2 P{w 2 )e~ sw \ (11) 
Jo 



which by (0) gives 

N 



(12) 



The cumulant generating function ^(s) can be obtained 
by expanding the logarithm as 



k V crT 1 -" / ' n u 
k=l v 7 n=l 



k N 



*(*)=lnG( S )=W(--y ( 13 ) 



Hence the fc-th cumulant of the roughness is 

N 

(14) 



n=l 



Note that for ak > 1 the cumulants can be expressed in 
terms of Riemann's zeta function in the limit of large N 
as 



(fc-1)! 
(o-T 1 -") 



-C(ak). 



(15) 



The low-order cumulants have the usual meaning of av- 
erage and variance 



Ki = (w 2: 



n— 1 



(16a) 



JV 



n 2 = (w 2 2 )-{w 2 ) 2 = 2^^- ( 16b ) 

The PDF corresponding to the generator ( [l2] ) will be 
the main concern of this paper. We shall also discuss 
"window" boundary condition, WBC, so in case of am- 
biguity the above quantities, pertaining to PBC, will get 
a subscript like G p (s). 



C. Scaling of averages 

The generating function of the PDF of the rough- 
ness still contains two scales, the observation time T 
(multiplied by cr 1 /*- 1 " - 1 ) and the microscopic time unit 
t = T/N. If N diverges, physical quantities will ex- 
hibit various scaling for different as. Whereas the above 
model is a very simple one from the viewpoint of the sta- 
tistical mechanics of critical phenomena, as it is Gaus- 
sian and massless, it is worth summarizing the scaling 
properties of averaged quantities, because they contain 
important information about the PDF. 

The fc-th roots of the fc-th cumulants ( |l4| ) for all ks 
have the dimensionality of the roughness, however, their 
scaling in the length T of the time interval may be dif- 
ferent. For k = 1, 2 these are the average and the mean 
deviation therefrom, respectively. Firstly, for a > 1 all 
cumulants converge in the limit N — > oo and we have 



H^k oc T Q ~ 



(17) 
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TABLE I: Scaling of the cumulants for various as. The dots 
indicate that the last formula is valid from then on, i.e., ". . . " 
extends the validity to higher ks and vertical dots to lower 

Qf-S. 



T"? 'A 71 O'O i~\T 

value of 01. 








(l,oo) 


rj-ia — 1 






1 


In ^ 

r 


0(1) 




(1/2,1) 


r a-l 


^IQ! — 1 




1/2 




V T 




(1/3,1/2) 




r a-l/2 jt-1/2 


yiti — 1 


1/3 






ji-2/3 3^ T T -2/3 


(1/4, 1/3) 






^_a — 1/3 rj-<— 2/3 rj-iCY — 1 



If a = 1 then the average, «i, logarithmically diverges in 
TV = T/t while all higher cumulants remain of unit order, 
independent of T and r. This corresponds to a critical 
dimension in second order phase transitions, but here we 
have a critical a = 1 where the logarithmic singularity 
occurs. For a < 1 the average no longer depends on T, 
rather it is proportional to r 1_Q , and for 1/2 < a < 1 the 
fc > 1 cumulants still exhibit the scaling ( |l7| ) . The second 
critical value of a is at 1/2, where the second cumulant 
has a logarithmic singularity. For a < 1/2 the be- 
comes proportional to j 1 - 1 / 2 T Q - 1 / 2 j while higher-than- 
-second order cumulants still satisfy ( Jl7| ) until the third 
critical value of a — 1/3 is reached. Continuing the rea- 
soning shows that at a^ = X/k the fc-th order cumulant 
develops a logarithmic singularity. A significant change 
occurs at a2 = 1/2, because for lesser as the scale of the 
mean deviation (fc = 2) will dominate over the scale from 
higher fc-s. This means that the PDF of the roughness 
on the scale of its mean deviation becomes Gaussian in 
the N — ► oo limit. Thereby the effect of the sequence of 
critical at-s for fc > 2 is suppressed. Table | summarizes 
the scaling of the cumulants. Note that the powers of T 
and r add up to a — 1 for each entry to produce the right 
time-dimensionality. 

As mentioned above and can be seen in Table ffl, for 
a < 1/2 the mean deviation goes with the — 1/2-th power 
of the size T. This result would follow if the central limit 
theorem could be applied to the roughness as given in Eq. 
(|J) . This is analogous to the mean field behavior beyond 
the upper critical dimension of a statistical mechanical 
system exhibiting a second order phase transition. 

Since the a — 1 value is a threshold in the sense that 
for a > 1, the scales of all cumulants diverge as T Q_1 , 
while for a < 1 the average becomes independent of T 
and higher cumulants vanish, the natural scaling for the 
two cases is different. If a > 1 the scaled quantity 



w 2 

Kl 



(18) 



has, in the limit N — > oo, a convergent PDF, devoid of 
any adjustable parameter. In the following, we shall refer 
to the use of variable x as scaling by the average. For 
a < 1 the same PDF would be a Dirac delta centered 
on 1, so one rather resorts to a scaling that effectively 
widens the delta peak. One can do that by introducing 

MM 



v 



W 2 - Kl 



(19) 



that is, scaling the roughness by the variance. The PDF 
as a function of y will clearly become Gaussian below 
a = 1/2, but develop a nontrivial shape for 1/2 < a. 
This variable was used in our recent Letter on the a — 1 
case Q, where we used the notation x for it. We will 
adopt the notation that functions of x and y carry the 
labels 1 and 2, respectively. So, e.g., the PDF of the 
scaled quantity ( |l8| ) will be denoted by 'E'i(x). 

The two scalings (|l8| ) and (|l9|) can lead to dramatically 
different shapes. For instance, for a — > 1 from above, Eq. 
( |l8| ) yields a peak centered at x = 1 and leads to a Dirac 
delta at x = 1 for a < 1. However, ( |l9| ) gives smooth 
functions continuously changing when a passes through 1 
and accordingly, we will use sometimes ( |i~9] ) in the region 
a > 1 for the sake of comparison. 



III. SCALING FUNCTION FOR PERIODIC 
BOUNDARY CONDITION 

A. Analytic approach for general a 

Knowing the explicit form of the generating function 
(|l2|), the roughness distribution for any finite N can be 
obtained by inverting the Laplace transformation ([ll] ) 



P p (a,w 2 ,T) 



ds 
2tt7 



N 

n 



(20) 



From here we concentrate on the N — > oo limit. In this 
limit the PDF of the scaling variable x or y, defined in 
Eqs. ( p^ , p^ ), becomes independent of T. Scaling by the 
variance, that is, using the y variable of Eq. (19), has the 
advantage to yield smooth PDF-s for all as. We refer 
to this type of scaling by the label 2 on the PDF. For 
a > 1/2 that PDF is 



$ 2p (a,y) = y/K2~P p (a, w 2 ,T) 



ds c y/C(2a)y S 

2iri 



Hi 



(21) 



The integrand has simple poles along the real axes and 
their contributions can be easily collected 



$ 2p (a,y) = VC(2qO 



m— 1 



m a e~ ma \ / ^ y ~ 1 Y 2 (a,m) , 

(22) 
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where 



Yi (a, to) 



n r 



n— l,7^m 



-(f)" 



(23) 



This series can be considered as the large-?/ expansion of 
the PDF. It is a general formula for any a > 1/2 and can 
be evaluated numerically. First, the precise numerical 
values of Y 2 (a,m) should be determined up to a certain 
w max . Greater m max is needed for smaller values of y. 
Since the sign of ^2(0;, to) is known to be (— l) m_1 , it suf- 
fices to evaluate \Y-2\. The logarithm of IY2I can be writ- 
ten as an infinite sum, which can be tackled numerically. 
Once the values of Y 2 (a,m) are given, the summation 
in (^2|) can be done easily due to the exponential cutoff 
in to. The resulting PDF is depicted on Figs. |l| and ^ 
for several values of a. On Fig. @ logarithmic scale was 
used for the ordinate to make the tails visible for several 
decades. It is clearly seen on these figures that the shape 
of the PDF depends smoothly on the value of a. 
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FIG. 1: Roughness distribution of periodic signals, Eq. (|22|), 
for different values of a as a function of the scaling variable 
([Til). Note that the PDF is a smooth function of a for all as. 

When a > 1, the natural choice for scaling variable 
is x of Eq. (filf), and the PDF in terms of x can be ob- 
tained by changing the variable in (|22] ) and making some 
simplifications 



$i p (a,x) = Ki$p(a,u7 2 ,7 1 ) 

OO 

= C(«) m a e~ xmaC(a) Yi(a,m) (24) 
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FIG. 2: The same as Fig. hi but on log-linear scale. 



to x. The PDF-s on Fig. [| were calculated this way. 
One can observe the singularity at x = 1 emerging as a 
approaches 1. 
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FIG. 3: Roughness distribution of periodic signals for differ- 
ent values of a as a function of the scaling variable ([l8]). Note 
that as a — > 1 the PDF dramatically sharpens and converges 
to a Dirac delta function centered at x — 1. 

For integer values of a > 1 the infinite product (|25|) 
can be written in a closed form 



a-l 



Y 1 {a,m)= A ' Y[ T(l - a k m) , (26) 

k=l 



a 



with 



OO 

Y 1 (a,m)= [] [l 



n— 1,^771 



(25) 



Note that it is more appropriate for numerical evalua- 
tion to use Eqs. (^2 23) and then change the variable y 



where a — exp(27ri/o). For a = 2 we recover the random 
walk result of Ref. Kal 



Yi(2,to) =2(-l)' 



(27) 



In the a = 4 case there are three factors in the product 
and using some basic properties of the T function one 
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obtains 



Fx (4, to) 



47rm(— 1 



,m — 1 



sinh(7rm) 



(28) 



a result which agrees with that of Ref . |33| . 

Using (p^), we can write the roughness distribution for 
any integer a > 1 in the relatively simple form 

&i P (a,x) = a((a) V ^ — i— ■ 

^— ' to! 

m— 1 

xe -m a c(«)x JJr(l-a fc m). (29) 
fe=l 

The curves for PBC in all figures with integer a > 1 were 
drawn based on this formula. 



Special cases 



1. 1// noise (a = 1) 



Here we briefly revisit the case of 1// noise presented 
in our recent Letter Jl | . The natural scaling goes now by 
the variance as in (|l9). From Eq. ( |2l| ) the generator of 
the scaling function $i p (l,x) is 



The formula (32b) is a special case of what is known 
as the Fisher-Tippet-Gumbcl function that emerges in 
extreme value statistics |28j. This comes about in 
a nutshell as follows. Suppose we have a random vari- 
able with some generic PDF, and we draw M times in- 
dependently from this distribution. The PDF of the A:-th 
largest of all those values, that is, the extreme value PDF, 
will be centered around a median increasing with M. 
Obviously, for M — > oc and k fixed, the extreme value 
PDF will be determined by the tail of the original PDF. 
In this limit, an appropriate linear rescaling of the fc-th 
largest value yields generally an M-independent scaling 
function for the extreme value PDF. Now, in practical 
terms, if the original PDF has a tail decaying faster than 
any power law, then the scaling function is the Fisher- 
Tippet-Gumbel distribution. Our scaling function fl32bj ) 
corresponds to the k = 1 case, scaled to have zero average 
and unit variance. 

The fact that ( |32b| ) is related to extremal value statis- 
tics does not reveal automatically the mechanism of selec- 
tion of extremes in the case of the roughness of a signal. 
Also, it should be emphasized that we do not have the 
usual extremal value distribution functions for as other 
than one. So our result raises the problems of why the 
1// noise is distinguished among all a-s and how extreme 
value selection comes about there rather than resolves 
them. 



G 2P (M)=n i 



,sf an 



(30) 



where the upper limit of the product could safely be taken 
to infinity and a = Tr/y6. This formula produces the 
Gamma function as 



G 2p (M) 



— P l s l 



° r ( 1+ ^) 

due~ u (ue 7 ) s/a . 



(31) 



where 7 is Euler's constant and we also displayed Euler's 
integral formula for the Gamma function |38| . Introduc- 
ing the variable x = — (lnu + 7) /a we finally get 



dxe- sx <S> 2p (l,x), 

-00 



<&2 P (l,x) = ae 



(32a) 
(32b) 



The inverse Laplace transformation on (|32a[) gives ( |32b[) , 
so $2p(l,a;) is the sought PDF of the rough ness scaled 
by the variance. Note that strictly speaking ( 32a ) is not 
a Laplace transformation anymore, due to the shift of 
the average of < &2p(l, %) to zero, but the inverse transfor- 
mation can still be performed. In such cases the Fourier 
transformation is better suited for the generating func- 
tion, because the variable y of Eq. (|l^) is not restricted 
to non-negative numbers, but we could still derive the 
scaling function (32b) within the Laplace formalism. 



2. Wiener process (a — 2) 

For the Wiener process with PBC the generating func- 
tion as well as the asymptotics of the PDF for small and 
large x has been derived in p2[ . Interestingly, the dis- 
tribution can be expressed in terms of a known function. 
The normalization (|l8|), natural for a > 1, will be used. 
For the integrated density with a = 2 we get from (E3) 



M lp (2,x) 



$i p (2,x) dx 



1 



2EK 

fc=i 



k e -r 



= 4 (O, e~* x/6 ). (33) 
i?4 is Jacobi's 4th Theta function fl38| a nd its second ar- 



gument is now the relevant variable 39| 

The PDF can also be written in a closed, implicit form 
in terms of complete elliptic integrals of the first and 
second kind (K(k) and E(k)): 

1 .(i 



<M2,*(fc)) 



x{k) 



k 2 f' A 

x K 3/2 (k) [K(k) - E(k)] (34 
6K(Vl ~ k 2 ) 



irK(k) 



(34b) 



The sum over poles described in section III A for a 



general a can be understood as the large- a; expansion 
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of the distribution. On the other hand, we also have 
a small- a; series for the special a — 2 case. This can 
be constructed by scaling the generating function Jl2] ) 
according to Eq. (|l8|), expanding it as 



IV. SCALING FUNCTION FOR WINDOW 
BOUNDARY CONDITION 

A. Simulations 



G lp (2,s) 



'6s 



sinh V6s 

oo 

k=a 



-(2fc+l)v^s 



and using the inverse Laplace transform 



c+zoo 



ds 
2ni 



a 2 - 2x 



-a 2 /Ax 



(35) 



(36) 



The result is 



*i P (2,aO = / — e-G lp (2, S ) 



,-3(2fc+l) 2 /2x 



(37) 



fe=0 



One recovers here the nonanalytic small- a; asymptotics 
shown in p2| , and one also has the corrections to it. 
Three terms from the sum suffice to produce the PDF 
<&ip(2,a;) with a uniform error bound of e = 4 x 10~ 5 in 
the sense that either the approximation is within e from 
the PDF, or the PDF is less than e. The precision of this 
approximation is, however, much better when the PDF 
is of order unity, to give a feeling of it one gets 29 dig- 
its of the PDF at x = 1. The series can be considered 
as the small- a; expansion of the derivative of the Thcta 
functi on (|33"1) , not availa ble in |38) . The formulas of Sec- 
tions [III B 1| and IIIB 2 can help to test the numerical 
evaluation of the series (E3j23). 



3. Large a limit 

For large a the lowest frequency mode dominates. In- 
deed, in the action (|J) the coefficients c n with n > 1 will 
have very small variance compared to the n — 1 case, 
so practically they are zero. The corresponding genera- 
tor function is obtained from (|12| ) by omitting all n > 1 
factors. Applying the scaling ( |18| ) we get 

G lp (oo,s) = (l+s)-\ (38) 

whence inverse Laplace transformation yields 



$i p (oo,x) = e 



(39) 



Recently in Ref. |35| a PDF was introduced and defined 
by cumulants for arbitrary dimensions. Its special case 
d = 1 corresponds to our cumulants in Eq. (]l4|). There 
the x 2 density with 2c? degrees of freedom for large a was 
found, that for d = 1 indeed gives (|39]) . 



Although periodic 1// Q signals exist |], [32| most 
of the experimental signals displaying 1// Q spectrum are 
not periodic. Having a long experimental signal, how- 
ever, the roughness can be calculated for small uncorre- 
cted segments (these are called the windows) and the 
roughness distribution can be constructed for this "win- 
dow" boundary condition, abbreviated as WBC. It is a 
plausible assumption that this PDF does not depend on 
the BC-s of the original signal provided that the size of 
the window, T, is much smaller than the size of the entire 
signal T p . Having uncorrelated windows would require 
large distances between them, however, in our experi- 
ence the PDF remained unchanged even for overlapping 
windows. In the simulations of this chapter always non 
overlapping, but neighboring windows were used. There- 
fore it appears that the PDF of a Gaussian, perfect 1//" 
signal with WBC can be computed having a long peri- 
odic signal of the same type. This method was applied 
numerically for general a and we report analytical results 
only for a = 2 and a — oo in Sec. [VB. 



The most accepted numerical way of generating a 
Gaussian l// a signal is to generate a Gaussian white 
noise first, filter the Fourier spectrum of it in order to 
get the desired l/f a behavior (i.e., after a fast Fourier 
transformation its real and imaginary parts are multi- 
plied by / _Q / 2 ) and finally transform it back. In this 
way the modes have amplitudes fluctuating around the 
desired l/f a value, and have random phases. Calculat- 
ing the roughness distribution of such periodic signals 
made it possible to check our theoretical predictions for 
periodic signals. The main advantage of generating pe- 
riodic 1//" signals is, however, that one can simulate 
WBC this way. 

Having the desired periodic signal of length T p one 
can construct easily the PDF of the roughness of non- 
overlapping parts of size T. The value of T p was chosen 
to be at least 2 20 while T varied between 2 6 and 2 18 . 
The PDF converged to a size independent shape for each 
values of a we have studied. The finite size effect was 
larger for smaller values of a. For a — 1 the PDF for 
T = 2 6 was already within a linewidth to the limit curve, 
however, for a = 1/2 such precision was only reached for 
T = 2 18 . 

The results for a few values of alpha are depicted on 
Figs. together with the two analytical results for 
WBC of Sec. IV B| , which provided a good test for the 
numerics. One observes that the curves are changing 
continuously with alpha in the whole range of [1/2, oo]. 
For a — 1/2 we recovered the Gaussian PBC result. It is 
also interesting to note that the a — > oo convergence is 
noticeably faster than that in the PBC case. For WBC 
already the a = 3 PDF can be well approximated by that 
for a = oo. 
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FIG. 4: Roughness distribution with WBC for different values 
of a as a function of the scaling variable (jl9|). Note that the 
a = 1/2, 2, and oo curves are exact results. 
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FIG. 6: Comparison of PBC (thick lines) and WBC (thin 
lines) for different values of a as a function of the scaling 
variable (p^|). 
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FIG. 5: The same as Fig. ^ but on log-linear scale. 

On Fig. H PDF-s with different BC-s are compared for 
several values of a. As we have already mentioned before, 
for a = 1/2 both BC result in a Gaussian PDF. For 
a = 1, as we already reported in Ref. Q, the difference 
between the PDF-s are relatively small, namely it was 
comparable to the precision of the experimental values. 
On Fig. I one can observe that the difference between 
the PDF with PBC and WBC are larger for larger values 
of a. 

We should emphasize that the thresholds a = 1 and 
a = 1/2 were found to play the same role for PBC and 
WBC. Namely, the natural scaling of the roughness goes 
by the average and by the variance for a > 1 and 1 > a, 
respectively, and for a < 1/2 the scaled PDF becomes a 
Gaussian for both BC-s. This evidence strongly supports 
our expectation that the scaling of cumulants for PBC, 
as summarized in Table |, also applies to WBC. 

It is worth recalling that critical exponents in sec- 
ond order phase transitions are generally believed not 



to depend on BC, but the scaling functions do vary be- 
low the upper critical dimension, the threshold for mean 
field behavior, for different BC-s @ §7], ||. The afore- 
mentioned BC-independence of the scaling of the cumu- 
lants on the one hand, and the BC-dependence of the 
scaled PDF for the roughness on the other one, show a 
close analogy with critical phenomena. The threshold for 
mean-field-type behavior is now a = 1/2, because below 
that the PDF is Gaussian. 



B. Analytic results in special cases 

1. Wiener process (a = 2) 

Since the increments of the Wiener process are uncor- 
rected, the trajectory in a window is coupled to the path 
outside only by its two endpoints. However, if the length 
of the outer trajectory is much larger than that of the 
window then the endpoints appear essentially as uncon- 
strained on the scale of the inner trajectory. Thus the 
WBC corresponds to the free BC in the a — 2 case. The 
generating function for free BC is, apart from scaling of 
s, the square root of the generating function for PBC as 
given below 



Giw(2,s) = 




V sinh \f\2s 

where the normalization Gi w '(2, s)| s= o = ~ 1 was used, 
corresponding to scaling by the average (|l^) . This yields 
the PDF <&i w (2,x) by inverse Laplace transformation. 
Eq. (^) can be understood if one notes that free BC 
means that the action is diagonalized by cosine eigen- 
functions, having zero derivatives at the endpoints. So 
the spectrum is labeled by integers, like for PBC, but 
there is no degeneracy. The generating function is in 
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essence the reciprocal square root determinant. While 
for PBC the square root disappears due to the two-fold 
degeneracy, in free BC the root remains. A more rigorous 
derivation will be presented in M] . 

The small- x series of the PDF $i w (2, x) can be ob- 
tained by expanding Eq. ( f40| ) as 

oo 

G lw (2, S ) = 2</T S V(-l) fc t 1 fe /2 )e-( 4fc+1 ^, (41) 



fc=0 

and using the Laplace transform of a single term in the 
sum as 



e-Hoo 



' S e 



ds 
2ni' 



a 3 / 2 



2 °V 4' 4' a 2 J 



,(42) 



where gFo is a hypergeometric function in standard no- 
tation |38[ . Based on these relations we finally get 



*iw(2,a:) 



c+zoo 



ds 
2tt7 



e sx G lw (2,s) 



x exp 



E( 4fc + 1 ) 7 2fcfc , 



fc=0 



3(4fc + l) : 



4.2: 

3 

"4 ; 



4j; 



3(4fc + l) 2 



(43) 



Although this can be considered as the small- a; expan- 
sion, it is enough to retain the first three terms to get 
an estimation of $i w (2,x) to a uniform precision of 
e = 6 x 10 -6 for any x. Similarly to the PBC case, this is 
understood such that either the approximation deviates 
from the PDF by less than e, or, the PDF is smaller than 
e. To illustrate the accuracy of three terms from (^) for 
reasonable x-s, at x = 1 they happen to give 53 digits of 
the PDF correctly. Series _fl43| ) has been used to draw the 
a = 2 curve of Figs. ||, ||, 



2. Large a limit 



As it has been discussed in Sec. [II B 3, for large a only 
the lowest frequency mode dominates the trajectory. In 
the full time interval with PBC, T p , this means that 



hit) 



2irt 



tp 



(44) 



where C\ = re lip /2 is the Fourier coefficient of the n = 1 
mode. The fact that the real and imaginary parts of C\ 
are independent, identically distributed, Gaussian vari- 
ables has the consequence that the PDF for the polar 
parameters is 



-r 2 /2a 



(45) 



where we do not specify the variance a as it will disappear 
from the final formula anyhow. 

Now we consider within the overall periodic signal h(t) 
a small window of length T -C T p . Since we shall average 
over the phase, we can take the window to be [0, T]. An 
expansion of ( J44| ) up to t 2 shows that the roughness of 
the trajectory within the window is in leading order 



W2 = r 2 6sin if, 



(46) 



where b is a constant depending on T and T p , whose 
value will turn out to be immaterial. We thus have for 
the PDF of the roughness 



pOC />27T 

P lw (oo,w 2 ) = dr I d<fp(r,if) 
Jo Jo 

xS (w 2 - r 2 6sin 2 tp) . (47) 

Inserting expression (|45| ) for p(r, tp) we get for the scaled 
variable x — w 2 / (w 2 ) the PDF 



$i w (oo,x) 



-x/2 



(48) 



Note that, similarly to the Wiener process, the large a 
limit of the generating function for the WBC is essentially 
the square root of that for PBC, apart from a scale change 
due to normalization. Indeed, the generating function for 
PBC is given by (p8|), while the Laplace transform for the 
above PDF is (1 + 2s)" 1 / 2 . 



l/f a SURFACES IN ARBITRARY 
DIMENSIONS 



So far we have considered random fields of one com- 
ponent that were functions of time, that is, a 1 + 1 di- 
mensional system. A natural generalization is to consider 
d + 1 dimensional random surfaces, where the substrate 
has d-dimensional coordinates x, and the field h[x) still 
has one component. Imposing PBC now means that the 
substrate is a d-dimensional torus. Again, we define the 
probability density functional of a surface through an ac- 
tion that depends on the Fourier components c n of the 



surface, n = (rii, . . . 114), rij 



-N, 



,N being integers, 



and c n = c*_ n . The spatial unit is £, the length of a 
period L — N£, the unit volume v — £ d , and the total 
rf-dimensional volume V — L d . (In the case of a usual 
surface d = 2 and V is the area of the substrate.) The 
action is 



S({c n }) = 2aV 1 -<*/ d J2'\n\ a 



(49) 



Here prime means that the summation excludes the ori- 
gin and counts only half of the remaining index space so 
that if a vector n is included then — n is not. For a = 2 
this is the action associated with the stationary distribu- 
tion of the Edwards- Wilkinson model [E2I. In the case 
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TABLE II: Scaling of the cumulants for various as for gen- 
eral dimension The dots indicate that the last formula is valid 
from then on, i.e., ". . . " extends the validity to higher k-s and 
vertical dots to higher d-s. 



Range or 










value of d. ki 


\f^2 








(0,a) v ald - 1 
a InE 

V 

(a, 2a) v"^" 1 


0(1) 

ya/d-l 








2a : 




y_ 

V 






{2a, 3a) 


v a/d-l/2 y- 


1/2 


ya/d-l 




3a 








y-2/3 



(3a, 4a) 



v a/d-\/3 y-2/3 ya/d-l 



of general cv-s, Eq. ( [49|) is the long-range interaction part 
of the single-component version of the generalized 0{n) 
Hamiltonian of jl3|, for a recent reference see 0]. We 
shall briefly review some scaling properties for arbitrary 
dimensions in order to put l// a noise, as discussed in 
previous Chapters, in a broader perspective. 

The roughness of a surface is a random variable, whose 
PDF c an b e derived similarly to the d = 1 case discussed 
in Ch. II B. The short-range interaction case, a = 2, has 



been studied in detail by Refs. ||24 
by (45) for arbitrary dimension. T 
of the PDF P(w2) is obtained as 



_ |3J for d = 2 and 

le generating function 



g( S ) = n 1 + 



aV i- a /d\ n \°< 



whence the cumulants are 



(fc-1)! 
{aV 1 - a / d ) 



\ak 



(50) 



(51) 



Note that the cumulants derived here were used to define 
the model of Ref. |jf|] . The sum converges for ak > d, 
diverges logarithmically as IniV oc la(V/v) for ak = d 
and like a power function as ]\[ d ~ ak oc (V/v) 1 ~ ak / d for 
ak < d. The scaling properties of the cumulants are 
summarized in Table |Dj. 

For fixed d, in the large a limit, the modes with \n\ = 1 
dominate and so the roughness W2 obeys the \ 2 distri- 
bution, as observed in Ref. pEfl . For finite a, the thresh- 
old dimension where the mean logarithmically diverges 
is d = a. For d > a, however, P(w2) becomes a Dirac 
delta, but if one looks at it on the scale of the variance for 
large but finite N then for d < 2a a nontrivial function 
emerges. At d = 2a the scale of the variance becomes 
larger than the scale of all higher cumulants, and thus for 
d > 2a the scaling function is Gaussian. This represents 



normal finite size scaling with fluctuations of the order 
of V^ 1 ! 2 . Table || is in accordance with the known fact 
that d = a and d — 2a can be viewed as the lower and 
upper critical dimension of the system, respectively [H . 

It is reasonable to assume that the scaling of the cumu- 
lants as described above is not specific to the PBC used 
here, but generally characterizes d+ 1 dimensional Gaus- 
sian surfaces with dispersion exponent a for any BC-s. 
Note that here only the scaling of averages were consid- 
ered, the evaluation of distribution functions in arbitrary 
dimensions is beyond the scope of the present paper. 



VI. FINAL REMARKS 

As we have shown, the roughness distribution of pe- 
riodic Gaussian 1/ f a signals can be calculated for arbi- 
trary a. The final expression is simple enough that it can 
be easily handled numerically and the scaling functions 
can be displayed in the relevant range of their argument. 
Also for WBC we provide a simple method to generate 
the scaling function by numerical simulation. Examining 
these scaling functions, we found an important feature in 
their a dependence. Namely, the shape of the functions 
varies noticeably with alpha in the physically rather in- 
teresting range of 1 < a < 2. This observation underlies 
the usefulness and effectiveness of the roughness distribu- 
tion as a tool for establishing common or distinct origins 
of scale- invariant behavior in different systems. 

The present gallery of scaling functions is ready to 
be applied for determining accurate values of a in scale 
invariant systems where the fluctuations are Gaussian. 
Since we have both the PBC and WBC scaling functions 
one can investigate models where PBC is used prefer- 
entially as well as experimental systems where WBC is 
usually obtained. Furthermore, the gallery can also be 
helpful in establishing the presence of non-gaussian ef- 
fects. It should be clear, however, the non-gaussian ef- 
fects are on the unfinished end of our study of roughness 
distributions. One can investigate the non-gaussian ef- 
fects in a given system by simulations pq , |3^] but the 
real question one should answer here is this: Can one 
find a classification of nonlinear theories which produce 
a given a, and can one find the roughness distributions 
for the various classes? Judging from the perspective of 
a related topic of critical dynamics this appears to be a 
highly nontrivial question. 
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